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Abstract 

The  stability  of  the  flow  of  a  viscous  incompressible  fluid  over  a  curved  compliant  wall  to 
longitudinal  Gortler  vortices  is  investigated.  The  compliant  wall  is  modeled  by  a  particularly 
simple  equation  relating  the  induced  wall  displacement  to  the  pressure  in  the  overlying  fluid. 
Attention  is  restricted  to  the  large  Gortler  number  regime;  this  regime  being  appropriate  to 
the  most  unstable  Gortler  mode.  The  effect  of  wall  compliance  on  this  most  unstable  mode 
is  investigated. 
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I.  INTRODUCTION 


Our  concern  is  with  the  effect  of  wall  compliance  upon  the  stability  characteristics  of 
the  flow  of  an  incompressible  fluid  over  a  wall  of  variable  curvature.  Recent  theoretical 
work  on  the  effects  of  wall  compliance  on  the  stability  of  fluid  flows  has  been  restricted  to 
flows  over  flat  compliant  walls  (here  by  flat  we  are  referring  to  the  unperturbed  state  of 
the  wall).  This  work  has  incorporated  various  models  of  the  resulting  wall  displacements 
eg.  viscoelastic  walls1  and  spring  backed  walls2,3  (both  isotropic  and  anisotropic).  Conse¬ 
quently,  this  work  has  considered  various  aspects  of  the  stability  of  Tollmien-Schlichting 
waves  and  the  efFect  that  the  underlying  wall  motion  has  on  their  stability  characteristics. 
A  review  of  these  results  can  be  found  in  Carpenter4. 

The  effect  of  wall  compliance  on  flows  over  curved  walls  has  recently  been  considered 
experimentally  by  Yurchenko,  Babenko  &  Koslov5.  They  restricted  their  attention  to 
centrifugal  instabilities,  namely  Gortler  vortices.  Although  no  precise  measurements  of  the 
growth  rates  of  the  initial  disturbances  have  been  presented  by  these  authors  their  results 
for  a  compliant  wall,  in  which  they  plot  regions  of  the  Gortler  number/ wave  number  space 
corresponding  to  maximum  amplification  rates  of  the  initial  disturbances,  show  a  shift  to 
the  left  as  compared  to  the  results  obtainable  for  the  a  rigid  wall.  These  results  then 
suggest  that  the  effect  of  the  compliant  wall  has  a  stabilizing  effect  on  the  growth  rate  of 
the  longitudinal  vortex  structure. 

In  the  case  of  a  rigid  wall  it  has  been  shown  by  Denier,  Hall  &  Seddougui6  that 
there  is  a  most  unstable  Gortler  mode  (ie  a  mode  with  maximum  growth  rate)  in  the 
parameter  regime  G  1  (here  G  is  the  Gortler  number)  and  has  wavenumber  given  by 
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k  =  0(G1/5)  with  growth  rate  0(G3f 5).  Such  a  result  was  found  by  considering  both  the 
inviscid  Gortler  modes,  found  in  the  parameter  regime  G  »  1  with  wavenumber  k  =  0(1), 
and  the  right  hand  branch  Gortler  modes,  found  in  the  parameter  regime  G  1  and 
wavenumber  k  =  0(G1^4).  By  their  very  nature  the  inviscid  modes  do  not  feel  the  effect 
of  the  wall  and  so  we  expect  that  these  modes  are  virtually  unaltered  by  the  presence  of 
a  compliant  wall.  In  the  case  of  the  right  hand  branch  modes  it  is  known  from  the  work 
of  Hall7  that  these  modes  are  confined  to  a  thin  region  internal  to  the  flow  and  away  from 
the  wall.  Hence,  we  expect  that  these  modes  are  also  unaffected  by  the  presence  of  the 
compliant  wall.  However,  in  the  case  of  the  most  unstable  Gortler  mode  it  is  known6  that 
these  modes  are  confined  to  a  thin  wall  layer  of  thickness  0(G~1/5).  We  anticipate  that 
these  modes  will  be  strongly  affected  by  the  presence  of  the  compliant  wall.  For  these 
reasons  we  will  confine  our  attention  to  the  most  unstable  Gortler  modes. 

The  outline  of  the  paper  is  as  follows.  In  §2  we  formulate  the  problem  under  consid¬ 
eration  and  derive  the  governing  equations  for  the  fluid  motion.  The  equation  of  motion 
for  the  wall  displacement  is  assumed  to  have  a  particularly  simple  form  (thus  facilitating 
the  analysis  to  be  presented  in  §3);  further  refinements  to  the  governing  equation  for  the 
wall  displacement  can  be  made  to  take  into  account  the  effect  of  fluid  substrates,  elasticity 
of  the  compliant  wall  etc4.  Such  issues  will  not  be  pursued  in  this  paper.  Subsequently 
in  §3  we  will  restrict  our  attention  to  the  high  Gortler  number  regime  and  consider  the 
wavenumber  regime  k  =  0(G ^5).  The  eigenvalue  problem  for  the  most  unstable  Gortler 
mode  will  be  derived  and  in  §4  we  present  some  results  obtained  from  a  numerical  solution 
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of  this  problem.  The  implications  of  these  results  to  the  stability  of  flow  over  a  curved 
compliant  wall  will  be  discussed. 


II.  FORMULATION  OF  THE  PROBLEM 

We  consider  the  steady  flow  of  a  viscous  incompressible  fluid  over  a  compliant  wall 
of  variable  curvature.  Assuming  that  L  is  a  typical  lengthscale  over  which  the  curvature 
changes,  U ^  is  a  typical  flow  velocity  a  great  distance  from  the  wall,  and  v  is  the  kinematic 
viscosity,  we  define  a  Reynolds  number,  Re,  by 


(1) 


and  suppose  that  with  respect  to  Cartesian  axes  x* ,  y* ,  z *  the  unperturbed  position  of 
the  wall  is  defined  by 

y*  =  LRe-1/2g(x'/L).  (2) 

Upon  the  imposition  of  some  (suitably)  small  pressure  perturbation  the  position  of  the 
wall  is  now  defined  by 

y*  =  LRe~l/2  {y(x*/L)  +  6rj*  (x*  / L,  iZe1/2^*/L)}  ,  (3) 

where  Sy*  is  the  wall  displacement  whose  motion  is  assumed  to  be  governed  by8 

TRe-1/2<550vV  =  ~6p\-  (4) 

Here  B o  is  the  flexural  rigidity  of  the  compliant  surface  (we  have  assumed  that  the  com¬ 
pliant  surface  is  isotropic)  and  V4  is  the  usual  bi-harmonic  operator 

V4  =  d4.  +  2dl.  d?.  +  <94  . 

3 


The  term  6p\  is  the  small  pressure  perturbation. 

Defining  nondimensional  variables 

{x,y,z)  =  (X*  ,Re1/2y\Re1/2z')/L, 

(u+,v+,w+)  =  (u*,Re1/2v',Rel/2w,)/U0O, 

p+  =  p' /(pU2  ),  p'  =  Pq  +  6p\, 
we  restrict  attention  to  the  limit  Re  — >  oo  and  write 

v"*",  w+)  =  (a,  a, 0)  +  S(u,  y,  w)  +  0(62), 

p+  =  p(x)  +  6Re~1/2p(x,y,z)  +  0(62),  pf  =  Re~1/2p. 

Here  the  basic  velocity  (a,®)  depends  only  on  x  and  y  whereas  the  velocity  perturbations 
( u,v,u> )  depend  on  all  three  dimensionless  coordinates.  Substitution  of  these  expressions 
into  the  Navier-Stokes  and  continuity  equations,  taking  the  limit  6  — ►  0  with  Re  held 
fixed,  equating  terms  of  order  6°,  and  taking  the  further  limit  Re  — >  oo  in  the  resulting 
equations  for  (a,C,0 ,p)  and  ( u,v,w,p )  gives,  to  leading  order  in  powers  of  .Re-1, 

a«  +  ®v  =  0,  aa.  +  ®ey  =  —pK  +  ayy,  (5) 


and 


a*  +  vy  +  w,  =  0, 


ait,  +  VUy  +  ?2ax  +  VUy  =  Uyy  +  U,  *  , 
ay.  +  VVy  +  UVm  +  VVy  =  —py  +  Vyy  +  y„, 

a  WZ  +  ®t&y  =  - P ,  +  Wyy  +  W,,  , 

together  with  the  governing  equation  for  the  wall  displacement 


(6) 


(7) 
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where  B  =  Re2Bo/(pU%>  L3)  is  a  nondimensional  constant  which  we  initially  take  to  be 

0(1). 


By  writing  down  the  Taylor  series  expansion  for  the  no-slip  conditions  about  y  =  g(x) 


we  obtain  the  boundary  conditions 

U  =  V  =  0  y  =  g(x), 

u  =  -nyf)(x,z),  v  =  -vyrj(x,z),  w  =  0  y  =  g(x). 
In  the  limit  y  — >  oo  we  have  the  requirement 


(8) 


u  —*  ue(x),  (u,v,w)^>0  y  — »  oo. 


(9) 


where  ut(x)  is  the  dimensionless  free  stream  velocity;  in  this  paper  we  will  confine  our 
attention  to  the  Blasius  profile  so  that  ut(x)  =  1.  In  order  to  highlight  the  effect  of  wall 
curvature  we  make  use  of  the  Prandtl  transformation 

y-*y  +  g,  v-tv  +  g'n,  v-*v  +  g'u, 

whilst  all  other  variables  remain  unchanged.  The  governing  equations  then  become 

Um  +  Vy  =  0,  GG*  +  GGy  =  —px  +  Gy  v  , 

(10) 

n  =  v  =  0  y  =  0,  G  — >  1  y  — >  oo, 

and 

w»  +  Vv  +  w,  =  0, 

UUX  4-  GUy  +  GO,  4-  VUy  —  uvv  4"  W,  ,  , 

Uvm  +  uGy  +  uVx  4-  Gvy  4-  G\{x)Uu  =  —pv  4-  vvv  4-  v««, 

UWS  +  VWy  =  ~ps  4-  Wyy  4"  W,  M  ,  (H) 

Bf)tlts  =  —  p, 

u  =  —Gv  77,  v  —  w  —  0  y  =  0, 

(u,  v,  w)  — ♦  0  y  —*■  00. 
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In  (11)  we  have  replaced  2 gxa  by  Gx(%)  where  G  and  \  will  be  refered  to  as  the  Gortler 
number  and  wall  curvature  respectively. 

Finally  taking  the  Fourier  transform  of  the  system  (11)  in  z,  eliminating  w,  p  and  the 


wall  displacement  77  from  the  resulting  equations  gives 


n 

Gu,  +  Gu„  +  u G,  +  vUv  =  uvy  -  k  u, 

{gsvv  +  k*  +  k2vv}  v  +  vxuyy  +  {uxxv  +  k2vx  +k2xGu}u 
4-  (  uvv  -  G-^y  +  k2u |  vx  4-  2  [  G*t.  +  G.  — }  ux 


vv  “dy2  '  n  “J  “  l  1  = 

4-  uyyyv  -  OiW  -  K  +  2*2}  vvv  +  {nxy  +  A:2g}  =  0, 

together  with  the  boundary  conditions 


(12a, b) 


-ft. 


u 


Bk6 


^  {—  (uxyv  +  Vyyy)  +  k2u*  +  u  ( a,  *  +  u* „ )  4-  v(uxy  4"  nvv)}  y  —  0, 


v  —  v  =  0  y  —  0, 


(13) 


u,v,v  — »  0  y  — >  oo. 


Here  fc  is  the  transform  variable,  (u,v)  are  the  transforms  of  (u,v).  For  0(1)  values  of  the 
wavenumber  fc  and  finite  values  of  the  Gortler  number  G  the  solution  of  the  system  (12,13) 


can  only  be  found  by  numerical  integration.  However,  in  the  asymptotic  regime  G  1 
progress  can  be  made  analytically  and  it  this  parameter  regime  which  we  consider  in  the 
following  section.  From  the  boundary  condition  (13)  we  see  that  the  case  of  the  rigid  wall 
is  recovered  in  the  limit  B  — ►  oo. 


Ill,  THE  most  unstable  gortler  mode 

In  the  case  of  a  rigid  wall  it  has  been  shown**  that  the  system  (12,13)  has  a  mode  with 
maximum  growth  rate  in  the  parameter  regime  G  >  1  and  the  wavenumber  of  this  most 
unstable  mode  is  k  =  0{G1/b).  Such  a  result  was  found  by  considering  the  structure  of 
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both  the  problem  of  the  inviscid  Gortler  modes,  which  are  found  in  the  regime  k  =  0(1), 
G  >  1,  and  that  of  the  right  hand  branch  modes,  found  in  the  parameter  regime  G  »  1 
and  k  =  0(G1/4). 

By  inspection  of  the  structure  of  both  the  inviscid  modes  and  the  right  hand  branch 
modes9  it  is  clearly  seen  that  in  the  case  now  under  consideration  these  modes  are  essen¬ 
tially  unchanged  from  those  relevant  to  the  rigid  wall  problem.  In  fact,  an  analysis  of  the 
system  (12,13)  shows  that  in  these  two  parameter  regimes  the  effect  of  wall  compliance 
is  encountered  as  a  lower  order  effect  and  does  not  contribute,  at  leading  order,  to  the 
structure  elucidated  in  Ref.  6.  This  result  is  easily  seen  if  one  considers  the  position  of 
both  the  inviscid  and  right  hand  branch  modes;  the  inviscid  modes  (by  definition)  do  not 
feel  the  effect  of  the  wall  whereas  the  right  hand  branch  modes  are  situated  at  a  position 
internal  to  the  flow  regime  and  so  will  not  feel  the  effect  of  the  compliant  wall. 

For  these  reasons  we  will  restrict  our  attention  to  effect  of  wall  compliance  on  the 
most  unstable  Gortler  mode.  To  consider  this  mode,  we  note6  that  this  mode  is  confined 
to  a  wall  layer  of  thickness  0(fc-1)  where  the  wavenumber  scales  as  k  =  G(G^5).  We  then 
write 

<p  =  ky,  k  =  A  G1/5.  (14) 

We  expand  the  mean  flow  in  this  layer  as 

H(x)k  V  +  ---> 
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and  expand  the  velocity  field  ( u,v )  and  the  surface  displacement  r?  as 


«  =  {lTo  +  G-1'5Ul  +  •••}  exp 
v  =G2/5{v0  +  G~1/5\\  +  •  ■•} 


/ 


.L 


(15) 


where  we  have  anticipated  growth  rates  of  0(G3^5)  (see  Ref.  6).  Substitution  of  (14)  and 
(15)  into  the  governing  equations  (12a, b)  we  obtain,  to  leading  order, 


{ 


f2 

dp 

dip'2 


X3 

A3 


XHpUo 
A3  • 


(16) 


An  evaluation  of  the  boundary  conditions  (13)  shows  that  the  effect  of  wall  compliance 
first  becomes  important  when  the  compliance  parameter  B  =  G(G~1/5).  Hence,  we  define 


Bo  =  BG1/5. 


The  boundary  conditions  appropriate  to  (16)  are  then  found  to  be 


u  d3V0 

Uo  =  J^7^’  v°  = °-  ’r"-°  ^  =  0’ 

Uo  =  0,  Vo  =  0,  Vq  =  0  <p  =  oo. 


(17) 


As  noted  above,  the  case  of  the  rigid  wall  is  recovered  in  the  limit  Bo  —*  oo. 

For  computational  purposes  it  is  convenient  to  eliminate  the  parameters  n,  \  from 
(16)  by  writing 


A  =  (x/i2)1'5*,  Bi  =  BoUv2)11*,  U0  =  »v 0.  Vo  = 
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in  which  case  (16,17)  may  be  rewritten  as 


{  d 2  /?*>  i  -  _  Vo 

t  dp2  A3  1/£/0_A2’ 


/V  ,1 

{  d 2  l) 

l  dp 2 

A3  / 

1  dp*  li 

(18) 


Uo  = 


dr  Vo 


\3B\  dp3  ' 

Uo  =  Vo  =  V'o  -  0 


V  0  =  Vn  =  0 


p  =  oo. 


P  =  0, 


For  a  given  value  of  the  parameter  B\  the  system  (18)  constitutes  an  eigenvalue  problem 
for  the  growth  rate  (3  =  /3(A). 

Before  proceeding  with  a  discussion  of  the  results  of  our  numerical  integration  of  the 
system  (18)  we  note  that  the  limiting  inviscid  form  of  the  governing  equations  is  recovered 
by  taking  the  limit  A  — >  0  while  the  limiting  form  of  the  right  hand  branch  governing 
equations  are  obtained  by  taking  the  limit  A  — ►  oo.  In  fact,  from  (18)  it  is  easily  seen  that 


/3~A1/2,  A  — +  0,  /3  ~  A-2,  A  — »  oo, 


thus  demonstrating  that  for  some  intermediate  value  of  A  there  exists  a  most  amplified 
mode. 


IV.  RESULTS  AND  DISCUSSION 

The  system  (18)  was  solved  numerically  using  a  finite  difference  scheme  in  the  in¬ 
dependent  variable  p.  For  a  given  value  of  the  parameter  B i  and  a  fixed  value  of  the 
wavenumber  A  an  initial  guess  for  the  growth  rate  0  was  supplied  and  the  program  was 
run  until  convergence  was  achieved.  The  eigenvalue  problem  (18)  was  then  solved  by 
continuously  iterating  in  the  wavenumber  A  to  generate  eigenvalues  (3  =  /3(A;i?i).  This 
procedure  was  repeated  for  various  values  of  the  parameter  B\. 
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For  a  given  value  of  B\  the  system  (18)  has  an  infinite  number  of  eigenvalues  (5n  = 
(3n(\;Bi).  In  Figure  1  we  present  a  plot  of  the  first  eigenvalue  (3  as  a  function  of  A  for 
various  values  of  the  parameter  B\\  these  eigenvalues  correspond  to  the  most  unstable 
Gortler  mode.  In  Table  I  the  maximum  eigenvalue  /3max  is  given  together  with  the  corre¬ 
sponding  critical  wavenumber  Acr.  We  see  that  the  introduction  of  wall  compliance  has  a 
stabilizing  effect  in  that  the  maximum  value  of  the  growth  rate  is  reduced.  From  Figure  1 
we  see  that  as  the  wall  compliance  is  increased  (ie  B\  is  decreased)  the  maximum  growth 
rate  decreases  until  a  threshold  value  is  reached.  Below  this  threshold  value  of  B\  the 
stabilization  due  to  the  wall  compliance  remains  constant;  from  Table  I  we  see  that  the 
stabilization  is  5%  compared  to  the  case  of  the  rigid  wall. 

From  Table  I  we  see  that  as  the  compliance  parameter  is  decreased  from  the  rigid  wall 
value  there  is  a  slight  increase  in  the  critical  wavenumber  until  a  threshold  value  of  the 
parameter  B\  is  reached  after  which  point  the  critical  wavenumber  decreases  in  magnitude. 

In  figures  2  and  3  we  present  the  eigenfunctions  for  the  maximum  eigenvalues  presented 
in  Table  I.  The  effect  of  wall  compliance  is  to  decrease  the  amplitude  of  the  eigenfunctions. 
Otherwise  there  is  little  qualitative  change  from  the  case  of  the  rigid  wall. 

In  conclusion,  the  presence  of  a  compliant  wall  has  the  effect  of  stabilizing  the  longi¬ 
tudinal  vortex  motion  which  is  in  qualitative  agreement  with  the  available  experimental 
results5.  However,  this  stabilization  is  relatively  small  and  hence  we  may  conclude  that 
the  stability  of  the  flow  over  a  curved  compliant  wall  is  essentially  unchanged  from  that  of 
the  case  of  a  rigid  wall.  It  remains  an  open  question  as  to  whether  either  finite  amplitude 
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effects  or  the  presence  of  wall  modes  (ie  modes  induced  by  the  motion  of  the  compliant 
wall)  will  significantly  influence  the  structure  discussed  in  this  paper. 
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